Total cross section for p — d breakup below 30 MeV 
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Abstract 

The total cross section for p-d breakup is studied in terms of the elastic 
iS-matrix through the unitary condition. Calculations using the complex 
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. Kohn variational method along with the Pair Correlated Hyperspherical Har- 

|Q . monic basis are presented. The results have been restricted to energies below 

E p = 30 MeV where Coulomb effects are expected to be sizable and are com- 
pared to the existing data. Two different measurements have been found in 

. 

the literature: 40 years ago, Gibbons and Macklin (1959); and 26 years ago, 
Carlson et al. (1973). The calculations are found to be in reasonable agree- 
ment with these old data, though a discrepancy is observed near the deuteron 
breakup threshold. Moreover, a detailed analysis of the contributions to the 
observable from different partial waves has been performed. Unexpectedly, 
the main contribution for a wide range of energies has been detected in the 
J = 3/2~ state. 
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Studies in the three-nucleon (3N) continuum are based mainly in N-d scattering experi- 
ments. The deuteron is the only existing bound state in the two-nucleon (2N) system with a 
binding energy Bd = 2.225 MeV. Therefore, for incident nucleon energies below 3.337 MeV, 
or incident deuteron energies below 6.675 MeV, there are no open channels (disregarding the 
very low probability for radiative capture or bremsstrahlung) and the reaction goes through 
the elastic channel. Accordingly the elastic scattering matrix, the iS-matrix, is unitary. 

For energies above the deuteron breakup threshold (DBT) three free nucleons are present 
in the outgoing channel. For example, in p-d scattering a free neutron can be observed only 
above the DBT (neutron production) whereas for the n-d reaction the observation of a free 
proton is only possible above the DBT. The elastic (S-matrix is no longer unitary and the 
missing flux in the elastic channel is related to the total breakup cross section. 

The total breakup cross section accounts for all possible configurations of the three free 
outgoing particles and its expression is given for example in Ref. |IJ in terms of the T-matrix 
elements connecting the elastic channels to the inelastic ones: 

^-^F (2/ 1+1 )W + i) ^ (2J+1) 5^ r (1) 

The index a labels the elastic channels while f3 runs over all possible configurations of the 
three outgoing nucleons. In addition k 2 = fr^o, where \x is the reduced nucleon-deuteron 
mass, Eq is the center of mass energy, and I\ (I2) is the spin of the incident particle (target). 
For each state with total angular momentum J the flux conservation in channel a imposes 
the unitary condition 

EI^I 3 + EW = i- (2) 

a' $ 

Experimental measurements of the total breakup cross section 07, are scarce and have 
been performed many years ago with limited accuracy. For the n-d system an indirect 
determination is possible through the subtraction of the total elastic cross section a e \{n — d) 
from the total nuclear cross section a to t{n — d). Existing a e i(n — d) data lack the desired 
accuracy producing a large uncertainty in cr&(n — d), particularly at low energies 0. Direct 



measurements of <Jb{n — d) have also been done more than 20 years ago ||. The interest at 
that time was to produce data for comparison with the first theoretical attempts to solve 
the 3N continuum using semi-realistic N-N potentials. These early potential models gave a 
qualitative description of the data but the uncertainties of ~ 10% were too large to make 
definitive conclusions. The present status of this observable is briefly given in Ref. || (see 
page 184) where modern theoretical predictions are compared to the data. The calculations 
in Ref. [f§] have been performed by solving the Faddeev equations in momentum space using 
different realistic potential models. The discrepancies observed at low energies have been 
attributed to the low quality of the data. In fact, since the total and elastic cross sections 
are well described by theory there is no reason to expect a discrepancy in the total breakup 
cross section. On the other hand, since the the total, elastic, and breakup cross section are 
related by the unitary condition and the optical theorem, it is highly desirable to have a 
consistent description of all three observables. New determinations of cr&(n — d) could come 
from both direct or indirect measurements. 

The situation of the p-d total breakup cross section is somewhat different. In this case 
no indirect measurement is possible due to the Coulomb divergence in the elastic amplitude. 
Two direct measurements have been found in the literature. A very old one by Gibbons 
and Macklin (1959) below E p = 5.6 MeV ||, with systematic uncertainties estimated to 
be 40% at E p = 3.5 MeV and decreasing to 5% at E p = 5.5 MeV. A second experiment 
was performed by Carlson et al. (1973) for proton energies between 22 and 46 MeV ||, 
with uncertainties between 8% and 5%. The experimental techniques utilized were quite 
different in the two measurements. In the first aj, has been determined by measuring the 
total neutron yield while in the second one an attenuation technique has been used. These 
two measurements provide valuable information to be used in the comparison to a theoretical 
description of p-d scattering including the long-range Coulomb force. In particular Coulomb 
effects are expected to be important in the Gibbons and Macklin experiment performed at 
energies near threshold. For example o"b(n — d) = (34 ± 6) mb at E n = 4.9 MeV (Holmberg 
(1969) [§]) to be compared to the measurement of Gibbons and Macklin, a^p—d) = 21.5±1.4 
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at the same energy. The Coulomb force reduces the breakup cross section by approximately 
30%. At higher energies the effect is less evident: at nucleon energies of 22 MeV the error 
bars of the n-d and p-d measurements start to overlap, and the effect seems to be extremely 
small above 40 MeV. A collection of n-d and p-d data is given in Fig. 1. 

Recently the variational technique based on an expansion in terms of the Pair Correlated 
Hyperspherical Harmonic (PHH) basis has been extended to describe elastic N-d scattering 
above the DBT [|7|-f| . The method consists of an expansion of the 3N wave function onto the 
PHH basis and the elastic elements of the iS-matrix are obtained through the application 



of the Kohn variational principle in complex form (see Ref. |T(| and references therein). 
The inclusion of the Coulomb interaction in the context of the variational method using 
the PHH expansion gives no particular troubles provided the correct boundary conditions 
are imposed on the wave function in the asymptotic region. An extensive discussion of the 
asymptotic form of the wave function in connection with the Kohn variational principle and 
the application of the PHH expansion is given in Ref. H]. In Ref. W\ it was shown that 



the method reproduces the benchmarks of Ref. [IT| for n-d scattering whereas in Ref. || 
the elastic cross section and polarization observables have been calculated for n-d and p-d 
scattering at nucleon energies of 5 and 10 MeV using realistic potentials. In the later case 
the applications to n-d scattering have been shown to be in close agreement to those given 
by the Bochum-Cracow group ||. 

Applications to describe p-d scattering above the DBT are of particular interest due 
to the historical difficulties to manage the distortion introduced by the Coulomb force in 
the asymptotic region. In this context, observables sensitive to the long-range Coulomb 
interaction give an unique opportunity to test techniques devised to take into account such 
distortion. The p-d total breakup cross section is a well suited observable to make compar- 
isons. The existing amount of data, though small, is sufficient to extract conclusion about 
the capability of the PHH technique to describe the charged 3N system above threshold. 
Moreover, as we shall see, it will be possible to evaluate the need for new and higher-precision 
measurements of <j\>. These measurements would increase our understanding of <Jb and im- 



pose stringent conditions on the theoretical methods and phase-shift analyses (PSA's). 

The breakup cross section c&(iV — d) is given in Eq. ([I]) in terms of the T-matrix elements 
connecting the elastic and inelastic channels. Using the fact that those elements are related 
to the elastic iS-matrix elements through the unitary condition of Eq. (0) it is possible to 
re- write the above equation in terms of the elastic matrix elements.. In the case of N-d 
scattering the elastic iS-matrix corresponding to a state with total angular momentum J 
is a 3 x 3 matrix. Each matrix element J S^'£' is labeled by the sets of quantum numbers 
[a = LS] coupled to J. Here L is the relative angular momentum between the incident 
nucleon and the deuteron and S is the total spin coming from the coupling of the spin of 
the deuteron Sj = 1 to the spin S^ = 1/2 of the nucleon. Accordingly S = 1/2 or 3/2 and 
there are three possible couplings of [LS] giving J and conserving parity [the parity of the 
state is given by (-1) £ ]. In the case of J = 1/2^ there are two possible couplings and the 
(S-matrix is a 2 x 2 matrix. Moreover for N-d scattering {21 \ + l)(2/2 + 1) = 6. From the 
above considerations Eq. ([]]) can be given in the more compact form: 

a b (N ~d) = ~ £(2 J + l)tr{Ij - SjS\} , (3) 

where Ij is the 3x3 identity matrix excepts for J = l/2 ± which is the 2x2 identity matrix. 
Eq.© gives 0\> in terms of the elastic matrix elements at fixed J, the sum runs over all 
possible values of J and parity (the sum over different parities is understood). Although 
the sum runs from J = to infinity there is a rapid convergence due to the fact that each 
Sj matrix becomes closer to unitary as J increase. High J values correspond to high L 
values and the centrifugal barrier prevents peripheral waves from participating actively in 
the breakup process. Moreover quartet states are also peripheral due to Pauli blocking. For 
each J states there are two quartet states and one doublet state which in combination with 
the increasing L values limit the number of terms in the sum to a tractable number inside the 
energy range we are going to consider. In fact, below 30 MeV states with J > ll/2 ± make 
a negligible contribution to Eq. (^) gives a very simple picture of the observable, as the 
diagonal elements of the product matrix SjS\ are always in the range (0, 1). Therefore the 
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quantity Aj = tr{Ij — SjSj} divided by tr{Ij} gives a measurement of the inelasticity in 
that state. The complete contribution to the breakup cross section is obtained after taking 
into account the spin degeneracy (2 J + 1). 

In the present paper we present calculations of <Jb{N — d) based on the PHH technique. 
Thorough details of the method are given in Refs. [0-0] and only a brief discussion will be 
presented. The starting point in the present calculations is the Kohn variational principle 
in complex form, which states that the functional 

i J S L L ' s S '} = J S L L ' S S ' + i(* Isj \H - E\n>s>j) , (4) 

is stationary with respect to the trial parameters in the three-nucleon scattering wave func- 
tion. For an incident state with relative angular momentum L, spin S and total angular 
momentum J it is 

^1sj= E [^(x i ,y i ) + ^ 5J (x i ,y i )] , (5) 

i=l,3 

and its complex conjugate is ^lsj- The ^ TS ^ term, ^q, describes the system when the 
three-nucleons are close to each other and, asymptotically, an outgoing three-particle state. 
Each amplitude ^/c( x j,y«), where Xj,yi are the Jacobi coordinates corresponding to the 
i-th permutation, has total angular momentum JJ Z and total isospin TT Z and it is decom- 
posed into channels labelled by angular-spin-isospin quantum numbers (/3-channels). The 
remaining two-dimensional amplitude is expanded in terms of the PHH basis 



x l1 y l )=p^ +L ^' 2 fp(x 



E4(p) (2) 4^( 



(6) 



L K 

where the hyperspherical variables are defined by the relations Xi = p cos <pi and yi = psin^j, 
fp(xi) is a pair correlation function and ^P^ L (4>) is a hyperspherical polynomial. The second 
term in the variational scattering w.f. describes the asymptotic motion of a deuteron relative 
to the third nucleon. It can be written in terms of the ingoing and outgoing solutions of the 
asymptotic N-d Schroedinger equation. 

^£Sj( X i> Yi) = ^£Sj( X i> Yd ~ E J ^LS S ^I/5'j( X i> Yi) ■ (?) 

L'S' 



For energies above the DBT the hyperradial functions u K (p) in Eq.(^) should describe an 
outgoing wave distorted by the Coulomb interaction between the two protons. Accordingly, 
for p — > oo, the hyperradial functions behaves as 

4(p)-E(^ xlos2Q 0?#'^e^ (8) 

/3>K' 

where Q 2 = ME/h 2 , with M the nucleon mass, E is the total energy (E = E — Bj), and 
the ^-matrix originates from the Coulomb potential. Extending the index (3 to include the 
hyperspherical quantum number K, the numerical constants B^ K are the T-matrix elements 
T J a fi defined in Eq.(|). 

We have studied the convergence of the quantity Aj in terms of the partial wave chan- 
nels used in the decomposition of the wave function for each state. Examples of the 
convergence can be found in Ref. ||. The accuracy of the results given here is estimated 
to be better than 1%. The calculations have been done with the Argonne vis (AV18) inter- 
action ]T2j and in some cases the three-nucleon force (3NF) of Urbana (UR) |13[ has been 
included for the study of the sensitivity to 3NF's. 

In Fig. 2 the PHH results for <Jb(p — d) are shown together with the experimental p-d 
data. The theoretical calculations (open squares) are in good agreement with the data. 
There is a window between 6 MeV and 22 MeV with no data, just in the region where the 
observable is varying rapidly and reaches its maximum. Above 20 MeV the calculations have 
been done at 22.7 MeV and 28 MeV, and are in close agreement with the measurements 
from Carlson et al. |J In Fig. 3 a detailed plot with semilogarithmic scale is given at the very 
low energies explored by Gibbons and Macklin ||. The theoretical results with the AV18 
interaction (open triangles) and AV18+UR interaction (open squares) are in agreement at 
the higher energies but start to deviate as the energy approaches the DBT. This region is of 
particular interest since Coulomb effects are huge and the breakup amplitude is attenuated 
by the Coulomb penetration factor exp(— 2tti]), with r] oc 1/^Eq — Bd- At 5 MeV the PHH 
n-d results are also given in Fig. 3 for the sake of comparison (diamond). The n-d data 
from Holmberg closest to 5 MeV are plotted also (open circles). The effect of the Coulomb 
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interaction is clearly evident. Moreover the calculations at 5 MeV based on the variational 
method reproduce both the n-d and p-d breakup reactions. 

Differences between theory and experiment are observed around 4 MeV and below. These 
discrepancies could originate from an incomplete treatment of the electromagnetic interac- 
tion. In fact, the calculations included the Coulomb interaction but not other small electro- 
magnetic parts of the interaction such as the vacuum polarization or the magnetic moment 
potentials. On the other hand systematic errors in the measurements are not unlikely - a 
new set of measurements will be highly desirable. 

In Fig. 3 the effects of the 3NF's are also studied (open squares). At these very low 
energies the breakup cross section is dominated by doublet states with L = 0, 1 since the 
probability to find the three nucleons close to each other is not reduced by either Pauli 
blocking or a high centrifugal barrier. The doublet 2 Si/2 state is sensitive to 3NF's. For 
example calculations with AV18 or AV18+UR at 3 MeV are 20% different for the 2 S\/2 



phase-shift In this case we are referring to energies below the DBT where the phase- 
shifts are real but we can expect sensitivity to 3NF's also in the imaginary parts above the 
DBT. The results given in Fig. 3 show a reduction of crj, of the order or 3%-4% when the 
AV18+UR interaction is used. The inclusion of the 3NF increases the binding energy of 
the three-nucleon system, and accordingly the inelasticity of the J = l/2 + state is reduced, 
diminishing the total breakup cross section. 

Finally we want to discuss the contribution of the different states to the total breakup 
cross section. Each state contributes to cr b an amount equal to (2 J + l)Aj. In Fig. 4 we 
plot this quantity up to J = 9/2 ± . It is interesting to note that above 6 MeV the state 
J = 3/2~ gives by far the main contribution to the observable. This state has the optimum 
combination of a large Aj value coming mainly from 2 P^/2 and the spin degeneracy 2 J+l = 4. 
This is twice the value for J = 1/2 and in fact the J = 3/2~ contribution is double the 
J = l/2~ one. Although the J = l/2 + state has the highest Aj value the degeneracy is too 
low and its contribution is the largest one only below 6 MeV. Above 24 MeV other states 
start to have important contributions, showing the weakness of the centrifugal barrier as 
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the energy increases. In Fig. 5 the contributions to cr& are given for energies below 7 MeV. 
The contributions calculated using the AV18+UR potential are also given (dotted line). As 
it was expected the J = l/2 + contribution is principal below 6 MeV and, more important, 
is the only one below 4 MeV. 

The reduced numbers of parameters in description of 0& at low energy is of particular 
interest. At 4 and 5 MeV a complete set of p-d vector and tensor analyzing powers as well 
as p-d differential cross section have been measured lqjl6 |. These high-precision data can 
be used to perform a single-energy PSA in which some of the phases must be allowed to 
be complex. From the above analysis it is clear that only few phases or mixing parameters 
have sizable imaginary parts. In addition, the inclusion of cr?, in the data base will impose 
further restrictions on the imaginary parts. 

PSA is extremely useful for performing comparisons between theory and experiment. 
Discrepancies at the level of observables can be traced back to single phases or mixing 
parameters which in turn are directly related to specific parts of the interaction. Above the 
DBT phases and mixing parameters are complex doubling the number of variables in the 
search procedure. Old PSA's encountered difficulties finding a precise determination of the 
real and imaginary parts of the phases and mixing parameters [17 . The high-precision p-d 
data obtained through the past years allowed for a better determination of the parameters 



through single energy PSA ||14]Jlq| . These analyses were limited to energies below the DBT. 
The study of the breakup cross section in terms of partial waves will be helpful when PSA's 
are extended to higher energies. In this light it would be particularly useful to perform 
new experimental determinations of 07, (p — d) at energies between 6 and 22 MeV, where 
there are presently no experimental results. It would also be most advantageous to perform 
these measurements at the same energies where highly-accurate differential cross section and 



analyzing power data already exist |TJj . 
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FIGURE CAPTIONS 

Fig. 1. Total breakup cross section for n-d (filled circles) and p-d (open circles) scattering 
up to 50 MeV. Experimental data are from Ref. []2],|3],|5||| . 

Fig. 2. Theoretical calculations of <Jb(p — d) up to 30 MeV (open squares). The experimental 
results of Gibbons and Macklin jS| (open triangles) and Carlson et al. || (open circles) are 
given for the sake of comparison. 

Fig. 3. Total breakup cross section below 6 MeV. Theoretical calculations for are given 
at four energies using the AV18 potential (p-d:open triangles, n-d:diamond) and AV18+UR 
potential (p-d:open squares). The experimental data are from Gibbons and Macklin || 
(p-d:filled circles) and Holmberg |3) (n-d:open circles). The solid line is a fit to the p-d data. 
Fig. 4. The quantity (2J+ l)Aj for different values of J and parity between 5 and 28 MeV. 
Calculations have been done at the same energies of Fig. 2. The solid lines are linear 
interpolations. 

Fig. 5. Same as Fig. 4 below 7 MeV. Calculations have been done using AV18 (solid line) 
and AV18+UR (dotted line). 
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